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ABSTRACT 

We present spectroscopic and photometric observations of the eclipsing sys¬ 
tem V1061Cyg, previously thought to be a member of the rare class of “cool 
Algols”. We show that it is instead a hierarchical triple system in which the 
inner eclipsing pair (with P = 2.35 days) is composed of main-sequence stars 
and is well detached, and the third star is also visible in the spectrum. We com¬ 
bine the radial velocities for the three stars, times of eclipse, and intermediate 
astrometric data from the HIPPARCOS mission (abscissae residuals) to establish 
the elements of the outer orbit, which is eccentric and has a period of 15.8 yr. 

We determine accurate values for the masses, radii, and effective temperatures of 
the binary components: MAa = 1.282 ± 0.015 Mq, i?Aa = 1.615 ± 0.017 R©, and 
= 6180 ± 100 K for the primary (star Aa), and M^h = 0.9315 ± 0.0068 Mq, 

Rxh = 0.974 ± 0.020 Rq, and = 5300 ± 150 K for the secondary (Ab). 

The masses and radii have relative errors of only 1-2%. Both stars are rotating 
rapidly {v sini values are 36 ±2 kms“^ and 20 ±3 kms“^) and have their rotation 
synchronized with the orbital motion. There are signs of activity including strong 
X-ray emission and possibly spots. The mass of the tertiary is determined to be 
Mb = 0.925 ± 0.036 Mq and its effective temperature is Tjl = 5670 ± 150 K. 
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The system is at a distance of 166.9 ± 5.6 pc. Current stellar evolution models 
that use a mixing length parameter auL appropriate for the Sun agree well with 
the properties of the primary, but show a very large discrepancy in the radius 
of the secondary, in the sense that the predicted values are ~10% smaller than 
observed (a ~5(T effect). In addition, the temperature is cooler than predicted 
by some 200 K. These discrepancies are quite remarkable given that the star is 
only 7% less massive than the Sun, the calibration point of all stellar models. 
Similar differences have been seen before for later-type stars, but the source of 
the problem has remained unclear. A comparison with the properties of other 
stars of similar mass as the secondary in V1061 Cyg has allowed us to identify the 
chromospheric activity as the likely cause of the effect. Inactive stars agree very 
well with the models, while active ones such as V1061 Cyg Ab appear systemat¬ 
ically too large and too cool. Theory provides an understanding of this in terms 
of the strong magnetic helds commonly associated with stellar activity, which 
tend to inhibit convective heat transport. The reduced convection explains why 
hts to models with a smaller mixing length parameter of ctml = 1-0 seem to give 
better agreement with the observations for V1061 Cyg Ab. 

Subject headings: binaries: close — binaries: spectroscopic — stars: evolution 
— stars: individual (V1061Cyg) —techniques: spectroscopic 


1. Introduction 

Accurately determined properties of stars in detached eclipsing binaries provide funda¬ 
mental data for testing models of stellar structure and stellar evolution (see, e.g., Andersen 
1991, 1997). For stars less massive than the Sun properties such as the stellar radius and 
the effective temperature have occasionally been found to disagree with model predictions 
(see, e.g., Lacy 1977; Popper 1997; Clausen et al. 1999a; Torres & Ribas 2002; Ribas 2003). 
Directed efforts to hnd additional systems in this regime suitable for testing theory (Pop¬ 
per 1996; Clausen et al. 1999b) have produced a few cases, while other examples have been 
found serendipitously (e.g., Creevey et al. 2005; Lopez-Morales & Ribas 2005). The present 
binary system is in the second category, since it was originally thought to be of a completely 
different nature. 

The photometric variability of V1061 Cygni (also known as HD 235444, HIP 104263, 
RXJ2107.3+5202, a = 21^07“2R52, 5 = +52°02'58M, J2000, SpT F9, V = 9.24) was 
discovered photographically by Strohmeier & Knigge (1959), and the object was classihed 
by Strohmeier, Knigge & Ott (1962) as an Algol-type binary with a period of 2.346656 days. 
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Other than occasional measurements of the time of primary eclipse, the system received very 
little attention until the spectroscopic work by Popper (1996), who observed it as part of his 
program to search for eclipsing binaries containing at least one lower main-sequence star (late 
F to K). On the basis of two high-resolution spectra and other information Popper concluded 
that V1061 Cyg was most likely a semi-detached system of the rare “cool Algol” class, and 
dropped it from his program. Unlike the classical Algols, which are composed of a cool giant 
or subgiant and an early-type star, the mass-gainer in the cool Algols is also of late spectral 
type (see Popper 1992). Since less than a dozen of these systems are known, V1061Cyg 
was placed on the observing list at the Harvard-Smithsonian Center for Astrophysics (CfA) 
in 1998 for spectroscopic monitoring, and photometric observations began later. Not only 
did we discover that it is not a cool Algol (it is well detached, as reported by Sheets, 
Marschall & Torres 2003), but we also found that: i) it is triple-lined (and a hierarchical 
triple); ii) the secondary in the eclipsing pair is less massive than the Sun and therefore 
potentially interesting for constraining models of stellar structure and evolution (Popper’s 
original motivation for observing it); and in) the mass ratio of the binary is quite different 
from unity, which makes it a favorable case for such tests. Furthermore, the comparison with 
theory shows a signihcant discrepancy in the radius of the secondary, corroborating similar 
evidence from other systems and providing some insight into the problem. 

We describe below our observations and complete analysis of this system, including a 
discussion of the possible nature of the deviations from the models for low-mass stars. 


2. Observations and rednctions 
2.1. Spectroscopy 

V1061 Cyg was observed at the CfA with an echelle spectrograph on the 1.5-m Wyeth 
reflector at the Oak Ridge Observatory (eastern Massachusetts), and occasionally also with 
a nearly identical instrument on the 1.5-m Tillinghast reflector at the F. L. Whipple Ob¬ 
servatory (Arizona). A single echelle order was recorded using intensihed photon-counting 
Reticon detectors, spanning about 45 A at a central wavelength of 5187 A, which includes 
the Mg I b triplet. The resolving power of these instruments is A/AA ~ 35,000, and the 
signal-to-noise (S/N) ratios achieved range from 13 to about 40 per resolution element of 
8.5 kms“^. A total of 74 spectra were collected from August 1998 until April 2005. One 
archival spectrum from 1987 obtained with the same instrumentation was used as well. That 
observation was taken by J. Andersen, possibly at the request of D. Popper (Andersen, priv. 
comm.), and it shows double lines (one sharp and one broad; see Figure 1). Popper (1996) 
referred to that observation, and remarked that his own two spectra of the star also showed 
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a sharp-lined and a broad-lined component. Most of our more recent observations display 
similar features. The fact that the sharp-lined star is hardly moving whereas the other one 
moves signihcantly led Popper to infer a mass ratio of roughly 8/1, and to conclude that 
VI061 Cyg is a cool Algol. 

As it turns out, the sharp-lined component is not the secondary in the eclipsing binary. 
Closer inspection of the cross-correlation functions, after we had estimated crude velocities 
and discovered that the sharp-lined and broad-lined stars were not moving together, revealed 
the presence of a third star (the true secondary) changing velocity with the same 2.35-day 
period as the broad-lined component. This is illustrated in Figure 1. In the following we 
refer to the broad-lined primary of the eclipsing binary (which is the more massive one of 
the pair) as star Aa, to the secondary as star Ab, and to the tertiary as star B. 

Radial velocities for the three stars were derived using an extension of the two-dimens¬ 
ional cross-correlation algorithm TODCOR (Zucker & Mazeh 1994) to three dimensions 
(Zucker, Torres & Mazeh 1995). This technique uses three templates, one for each star, and 
is well suited to our relatively low S/N spectra allowing velocities to be obtained reliably even 
when the spectral lines are blended. Templates for the cross correlations were selected from 
an extensive library of calculated spectra based on model atmospheres by R. L. Kurucz^, 
computed for us by Jon Morse (see also Nordstrom et ah 1994; Latham et al. 2002). These 
calculated spectra are available for a wide range of effective temperatures (Teg), projected ro¬ 
tational velocities (nsini), surface gravities (loggf), and metallicities. Experience has shown 
that radial velocities are largely insensitive to the surface gravity and metallicity adopted for 
the templates. Consequently, the optimum template for each star was determined from grids 
of cross-correlations over broad ranges in temperature and rotational velocity (since these are 
the parameters that affect the radial velocities the most), seeking to maximize the average 
correlation weighted by the strength of each exposure. Surface gravities of \ogg = 4.0 and 
logg = 4.5 were adopted for Aa and Ab, from the light curve analysis described in §2.3, 
and for star B we adopted loggf = 4.5. Solar metallicity was assumed throughout, which the 
discussion in §5 shows to be a good assumption in this case. The procedure of optimizing 
the templates leads to best-ht parameters for the stars as follows: 


Aa: 

TeS 

= 6160 ± 100 K ^ 

, nsini = 

36 ± 2 km s 

Ab: 

TeS 

= 5400 ± 200 K ^ 

, nsini = 

20 ±3 kms 

B: 

TeS 

= 5670 ± 150 K ^ 

, vsini = 

2±3 kms 


The stability of the zero-point of our velocity system was monitored by means of expo¬ 
sures of the dusk and dawn sky, and small run-to-run corrections were applied in the manner 


^Available at http: //cfaku5. cfa. harvard. edu. 
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described by Latham (1992). In addition to the radial velocities, we derived the light contri¬ 
butions of each star following Zucker, Torres & Mazeh (1995). These are f'Aa = 0.75 ± 0.01, 
^Ab = 0.09 ± 0.01, and £b = 0.16 ± 0.01, expressed as fractions of the total light, and they 
correspond to the mean wavelength of our observations (5187A). 

Due to the narrow wavelength coverage of the CfA spectra there is always the possi¬ 
bility of systematic errors in the velocities, resulting from lines of the stars moving in and 
out of the window with orbital phase (Latham et al. 1996). Occasionally these errors are 
signihcant, and experience has shown that this must be checked on a case-by-case basis 
(see, e.g., Torres et al. 1997, 2000). For this we performed numerical simulations in which 
we generated artihcial composite spectra by adding together synthetic spectra for the three 
components, with Doppler shifts appropriate for each actual time of observation, computed 
from a preliminary orbital solution. The light fractions adopted were those derived above. 
We then processed these simulated spectra with the three-dimensional version of TODCOR 
in the same manner as the real spectra, and compared the input and output velocities. The 
differences were typically small for stars Aa and B (under 1 kms“^), and were up to about 
3 kms“^ for star Ab. We applied these differences as corrections to the raw velocities, al¬ 
though they made little difference in the results other than slightly decreasing the residuals 
from the orbital hts. The hnal velocities including these corrections are given in Table 1. 
Similar corrections were derived for the light fractions, and are already accounted for in the 
values reported above. As a further test we changed the template parameters for the two 
fainter stars, which are more uncertain, and recomputed the radial velocities and corrections. 
The resulting velocities were hardly different, giving us conhdence in the robustness of the 
procedures. 

A double-lined spectroscopic orbital solution for the Aa-I-Ab system (see column 2 of 
Table 10 below) produced a reasonably good £t with rms residuals of 2.4 and 5.1 kms“^, at 
about the level expected based on the broad lines and faint secondary. However, the residuals 
showed an obvious downward trend, less obvious but still present also in the secondary. The 
velocities for star B show a drift in the opposite direction (Figure 2) with a slope about 
twice as steep, or ~5 kms“^ per 1000 days. The system is thus a hierarchical triple, with 
a period for the outer orbit that could be very long, judging from the data available. The 
spectroscopic material alone does not sufficiently constrain the wide orbit to properly account 
for the change in the center-of-mass velocity of the eclipsing pair, which could introduce biases 
in the derived masses of components Aa and Ab. However, as we show in the next section, 
the additional information provided by the measurement of times of eclipse and, to some 
extent, the astrometric data (ESA 1997) from the HIPPARCOS mission described later are 
very helpful in constraining the elements of the outer orbit. 
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2.2. Times of eclipse 

Measurements of the times of minimum light for V1061 Cyg spanning more than seven 
decades have been made by photographic, visual, and photoelectric/CCD techniques. Table 2 
collects all available estimates, including some that have not appeared previously in print 
and were kindly communicated to us by J. M. Kreiner (see Kreiner, Kim & Nha 2000). A 
total of 43 timings for the deeper primary minimum and 6 for the secondary are listed. Many 
of the more recent measurements between 2001 and 2004 are based on our own photometric 
observations with a robotic telescope described in more detail below. The O — C residuals 
from a linear ephemeris are shown in Figure 3, and exhibit a pattern of variations caused by 
the presence of the distant third component (light travel effect). This is most obvious in the 
more accurate and numerous measurements of the last 25 years, during which more than a 
full cycle has been covered. The period of the outer orbit is therefore fairly well constrained 
by these data, as is the amplitude of the variation, which is directly related to the velocity 
amplitude of the binary around the center-of-mass of the triple. These measurements provide 
information complementary to that from the radial velocities, and we incorporate them into 
the analysis in §3. 


2.3. Photometry 

Absolute photometry of VlOblCyg on the UBV system was obtained by Lacy (1992). 
The value of the reddening-free parameter Q = {U — B) — 0.72(5 — V) was calculated and 
compared with the standard star values of Johnson & Morgan (1953). The value corre¬ 
sponded to that of a spectral type GO star with an intrinsic color of (5 — V)o = O.OOlgiio, 
with no significant interstellar reddening. The uncertainty here comes mainly from the cal¬ 
ibration. The observed mean color index corresponds to a spectral type F9 star with a 
temperature of 6040 K, according to the calibration of Popper (1980). After some iterations 
with the photometric light curve htting algorithm described below, it was found that the 
visual surface brightness ratio JAh/Jha. is about 0.43, which constrains the difference in the 
visual flux parameter and hence the temperature difference through the calibration in 
Table 1 by Popper (1980). This procedure leads to photometric temperature estimates of 
6200 ± 100 K for the primary and 5280 ± 100 K for the secondary component. These tem¬ 
peratures were used to estimate the limb-darkening and gravity-brightening parameters to 
be used in the photometric modeling. They are in excellent agreement with those of the 
spectroscopic analysis above. 

Relative photometric observations of V1061 Cyg were obtained at two different facilities. 
Measurements in the V band were made with the 0.26-m URSA robotic telescope at Kimpel 



Observatory on the campus of the University of Arkansas at Fayetteville. A total of 6914 ob¬ 
servations were collected between March 2001 and October 2004, which are given in Table 3. 
For details on the observation and reduction methods the reader is referred to the description 
by Lacy, Claret & Sabby (2004). The precision of an individual measurement is estimated 
to be 0.010 mag. Observations were also made with the 0.4-m f/11 Ealing Cassegrain reflec¬ 
tor at Gettysburg College Observatory, in Gettysburg (Pennsylvania). The detector was a 
Photometries (Roper Scientihe) CH-350 thermoelectrically-cooled CCD camera with a back- 
illuminated SITe 003B 1024 x 1024 chip and standard Bessell BVRI hlters. Tables 4-7 list 
the observations obtained between September 2002 and September 2003, which number 732, 
738, 739, and 740 in the B, V, R, and / passbands, respectively. The estimated precision 
of these differential measurements is 0.011 mag in B, 0.008 mag in V, 0.009 mag in R, and 
0.012 mag in I. The comparison stars used for these two data sets are given in Table 8. 
The URSA photometry was referenced to the magnitude corresponding to the sum of the 
intensities of all four comparison stars in order to improve the measurement accuracy. 

Because the eclipse timings for V1061 Cyg described earlier show changes due to motion 
about the center of mass of the triple system, light elements used in analyzing the photometric 
data were based on contemporary observations made between March 2001 and October 2004 
with the URSA telescope. As seen in Figure 3 the curvature of the O — C residuals is 
relatively small in this interval, and a linear ephemeris is sufficient for our purposes. A £t 
to 12 primary minima (see Table 2) yielded 

Min I (HJD) = 2,452,015.90562(7) + 2.34663383(33) ■ E , 

where E is the number of cycles elapsed from the epoch of reference, and the values in 
parentheses represent the uncertainties of the elements in units of the last digit. Photometric 
data were then phased according to these elements. 

The URSA and Gettysburg data sets were analyzed with the Nelson-Davis-Etzel model 
(EBOP code) as described by Etzel (1981) and Popper & Etzel (1981). Although other 
models may be more sophisticated, this program is perfectly adequate for well-detached 
systems such as V1061Cyg. The main adjustable parameters in the NDE model are the 
relative surface brightness (J = JAh/Jxa) of the secondary of the eclipsing binary (star Ab) 
in units of that of the primary (Aa), the relative radius of the primary (rAa) in units of the 
separation, the ratio of radii {k = rAb/^^Aa), the inclination of the orbit of the binary (^a), the 
limb-darkening coefficients (xAa and XAb), the eccentricity parameters ecosci; and esinci;, and 
the amount of third light (^b, expressed as a fraction of the total light), which is significant in 
our case. The luminosity due to the reflection effect was computed from bolometric theory. 
The mass ratio (g = 0.7266) was adopted from the spectroscopic results in §3. Tests showed 
the orbital eccentricity was probably not significantly different from zero, consistent with 



the results described later, so e was fixed at zero for the photometric modeling. For the 
URSA data the EBOP algorithm converged with all variables free, except for the gravity- 
brightening coefficients y, which were set from theory (Alencar & Vaz 1997; Claret 1998). 
For the Gettysburg data sets the limb-darkening coefficients x had to be fixed at theoretical 
values (Al-Naimiy 1978; Wade & Rucinski 1985; Diaz-Cordoves, Claret & Gimenez 1995), as 
were the gravity-brightening coefficients. Additionally, for the Gettysburg data the values 
of third light in the spectral bands R, i?, and I had to be estimated from the temperature 
of the third star based on the spectroscopic analysis in §2.1. The third light value in V 
was adopted from the URSA data solution (£b = 0.144). For these reasons, the final mean 
geometric parameters based on the Gettysburg data have been assigned nncertainties two 
times the internal estimates, consistent with our experience with EBOP in these situations. 

Initial light curve solutions with the URSA data showed that observations during the 
night of JD 2,452,054 (24 May 2001), which followed a secondary eclipse, were significantly 
fainter than the fitted orbit by ~0.025 mag. No other observations at this phase (0.57) were 
collected until the 2003-2004 seasons. Additionally, data on JD 2,453,218 and JD 2,453,232 
(2004 season) obtained during and after secondary eclipse were seen to be brighter than the 
fitted curve by ~0.01 mag. Measurements taken 12 days earlier at a similar phase appear to 
be normal, as are those taken two months later. Possible explanations for these deviations 
include some unknown calibration problem in the data redactions, or transitory spots on 
one or both stars that change on timescales of a month or two. The latter would not be 
unusnal given the strong X-ray emission from the system (Voges et al. 1999), indicative of 
chromospheric activity (see §5.2). If due to spots, the extent of the active regions appears to 
be qnite large (up to 40°-50° in longitude), judging from the phase interval affected. Tests 
showed that removal of these URSA observations had an insignificant effect on the results 
other than decreasing the scatter. We have thus chosen to exclude the data on those three 
nights from the final fitting process. Purther snpport for intrinsic variations is given by 
similar deviations seen in the Gettysburg data at other phases, particularly phase ~0.45, as 
shown below in Figure 7. The angular extent of this particular feature is again ~40°. A 
hint of a dimming at abont the same phase may be present also in the URSA observations 
dnring the same season, althongh the data sampling is not optimal for confirming this. 

The results of our photometric analysis are given in Table 9, separately for each of the 
Gettysburg passbands and for URSA V. The mean of the Gettysburg results are also listed, 
with uncertainties that acconnt for the scatter between the BVRI results. There is excellent 
agreement in the geometric elements between the two data sets. The light curve and residnals 
for the URSA measurements are shown in Fignre 4, and enlargements of the regions aronnd 
the minima are shown in Figure 5 and Fignre 6. The Gettysbnrg light curves and residnals 
are shown in Figure 7. The binary is well detached: the stars occnpy only ~7% and ~2% 
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by volume of their critical equipotential surfaces. Both components are nearly spherical in 
shape (oblateness < 0.005, as dehned by Etzel 1981). The secondary eclipse is total, and the 
fraction of the primary light blocked at phase 0.0 is 40%. 

A common problem in solving the light curves of some eclipsing systems is the inde¬ 
terminacy in the ratio of the radii, particularly when the components are similar and the 
eclipses partial. As a test, we htted the URSA light curves with hxed values of k over a wide 
range around the best-fit value of 0.603. Solutions converged only between 0.55 and 0.63. 
Some parameters such as and to a lesser extent the relative brightness of the secondary 
(f'Ab), showed a very strong sensitivity to k. Only hts in the vicinity oi k = 0.603 (which 
minimizes the sum of residuals squared) resulted in values of £b and Ixh consistent with 
the spectroscopic estimates from §2.1 (after a small correction from the mean wavelength of 
5187 A to the V band). We take this as an indication that the adopted hts are realistic. 


2.4. Astrometry 

V1061 Cyg was observed by the HIPPARCOS satellite under the designation HIP 104263, 
from December 1989 through February 1993. The measured trigonometric parallax is tthip = 
6.25 ± 1.06 mas. The photometric variability of the star with the known 2.35-day period 
was clearly detected in the more than 120 brightness measurements made over the duration 
of the mission. Additionally, even though the wide orbit in V1061Cyg is much too small 
to be spatially resolved by HIPPARCOS, the astrometric solution revealed perturbations 
in the motion across the sky that led to its classihcation as a “Variability-Induced Mover” 
(VIM). These perturbations arise from the fact that the center of light of the triple system 
shifts during the eclipses due to the change in the relative brightness between the binary 
and the third star. The shift is approximately 25 milli-arc seconds, which the satellite was 
able to detect. As a hrst-order approximation the HIPPARCOS team assumed the angular 
separation between the binary and the third star was constant, so that the motion of the 
photocenter is along a straight line connecting them. With this they were able to infer a 
lower limit to the separation of approximately 70 mas, and a position angle (North through 
East) for the third star relative to the binary of 319° ± 14°. This linear and periodic motion 
of the photocenter was accounted for in deriving the position, proper motion, and parallax 
of V1061 Cyg, as reported in the catalog. 

In reality the third star and the binary are in orbit around each other, and from Figure 3 
the period appears to be of the order of 15 years. The 3.2-yr interval of the HIPPARCOS 
observations represents a non-negligible fraction of a cycle, and depending on the orienta¬ 
tion of the orbit this may have an impact on the derived parallax and proper motion, as 
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well as the inferred motion of the photocenter. In fact, turning the problem around, the 
HIPPARCOS observations can potentially provide valuable constraints on the outer orbit. 
Thus we incorporate these observations (in the form of the intermediate astrometric data, 
or “abscissae residuals”) into the orbital solution described next. 


3. Orbital solution 

Due to the much larger orbit of the third star compared to the eclipsing binary in 
V1061 Cyg, to first order we assume here that the hierarchical triple system may be separated 
into an inner orbit and an outer orbit, the latter treated as a “binary” composed of the third 
star (B) and the center of mass of the eclipsing pair (A). Orbital elements that refer to the 
outer orbit are indicated below with the subindex “AB”, and those pertaining to the inner 
orbit are distinguished with a subindex “A”. 

In addition to providing the spectroscopic elements of the inner orbit, in principle the 
radial velocity measurements yield also the elements of the outer orbit, except that in our 
case the coverage is insufficient for that purpose (see Figure 2), and the period of the outer 
orbit (Tab) is undetermined. Additionally, the semi-amplitudes of the velocity variation in 
the wide orbit (A'a and K-q) are poorly constrained because the observations do not cover the 
velocity extrema. The times of eclipse help in two important ways: they constrain the period 
(~15 years from a preliminary analysis), and they constrain the amplitude of the third body 
effect (Ko-c), which is directly related to Ka through Kq-c = A^aV"! — Cab/^ttc, where 
Cab is the eccentricity of the outer orbit and c is the speed of light. The fact that the 
inner binary is eclipsing allows the masses of the three stars to be determined. For this we 
incorporate the angle iA from §2.3 into the solution, along with its uncertainty. 

The use of the HIPPARCOS data introduces several other orbital elements into the 
problem that are constrained by the astrometry, including the inclination angle of the wide 
orbit (^ab) and the position angle of the ascending node (Dab, referred to the equinox of 
J2000). Additionally we must solve for corrections to the catalog values of the position of 
the barycenter (Aa*, AS) at the mean catalog epoch of 1991.25, corrections to the proper 
motions (A/i*, A/i^), and to the trigonometric parallax (Att)^. Because the inner binary 
is eclipsing there is redundancy in that the angle Ab can also be derived purely from the 
spectroscopic elements and the known value of A- We may thus eliminate one parameter, 
and for practical reasons we have chosen to eliminate and to retain Ab explicitly as an 
adjustable variable. 


^Following the practice in the HIPPARCOS catalog we define Aa* = AacosJ and A/x* = Afia cos <5. 



Altogether there are 19 unknowns in the least-squares problem: {Pa, Kas., K^h, ca, i^Aa, 
Ta} for the inner orbit, (Pab, A'a, ^ab, cab, ^a, f^AB, ^ab} for the outer orbit, the center 
of mass velocity 7 , and the HIPPARCOS-related elements (Ao;*, AS, A/i*, A/i^, Avr}. The 
longitudes of periastron u refer to the star or system indicated; Ta is the time of primary 
eclipse in the inner orbit, and Tab represents the time of periastron passage in the outer 
orbit. 

The amplitude of the motion of the photocenter in angular units, a'^^, is a function of 
the fractional mass and luminosity of the third star relative to the system as a whole. It 
is given by = a'A-Q^B — {3), in which the mass fraction is P = M-q/^Mas. + TfAb + M-q) 
and the luminosity fraction is /? = L-q/^Las. + Lxh + Tb). With the definition of third 
light used in EBOP, (3 = £b- The symbol represents the relative semimajor axis of 
the wide orbit expressed in angular measure. The notation a" is used here to distinguish 
from the semimajor axis in linear units, a. The mass fraction B can be expressed in terms 
of other elements, and Oab follows from knowledge of Pab, the masses, the parallax, and 
the application of Kepler’s Third Law. The luminosity fraction needs to be specihed in 
the HIPPARCOS passband (Pp), which is slightly different from Johnson V. To derive it 
we have made use of the HIPPARCOS epoch photometry for VlOblCyg, and solved the 
light curve with EBOP. Since the measurements available are fewer in number compared 
to our other light curves, we hxed the geometric elements to the average of the URSA and 
Gettysburg results from Table 9, and solved only for the passband-dependent elements J 
and £b- Limb-darkening coefficients have little effect and were hxed to their values for the V 
band. The ht to the HIPPARCOS photometry is shown in Figure 8 . The results we obtained 
from this adjustment are £b = 0.12 ± 0.02 and J = 0.44 ± 0.04 (Bp). With this value of 
£b the photocentric semimajor axis is completely determined at each iteration of the 
least-squares problem and is not an additional free parameter in the orbital solution. The 
formalism for incorporating the abscissae residuals from HIPPARCOS into the ht follows 
closely that described by van Leeuwen & Evans (1998) and Pourbaix & Jorissen (2000), 
including the correlations between measurements from the two independent data reduction 
consortia (ESA 1997). 

While the HIPPARCOS observations provide information that is complementary to that 
from other available measurements, the astrometric constraint on some of the orbital ele¬ 
ments is not particularly strong due to the limited precision of those measurements. Typical 
uncertainties for a single abscissa residual range from 2 to 4 mas, which is only a few times 
smaller than the motions we are trying to model. Additional constraints may be placed 
on the parallax by using the physical properties of stars Aa and Ab derived in §4, which 
allow the absolute luminosity of each component to be determined very accurately from their 
ehective temperature and radius. This is the same principle that has been used to great ad- 
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vantage for applications such as establishing the distance of external galaxies from eclipsing 
binaries (e.g., Guinan et ah 1998; Fitzpatrick et al. 2003). From the radii and temperatures 
in §4, along with the total apparent magnitude of V1061 Cyg and the fractional luminosities 
of the components, we obtain parallax estimates of 5.90 ±0.26 mas (Aa) and 5.72 ±0.48 mas 
(Ab), which are consistent with each other. We incorporate these estimates into the global 
solution as “measurements”, along with their uncertainties. 

The delay or advance in the times of eclipse caused by the third star is accounted for 
as described by Irwin (1952). The best-£t time of primary eclipse (Ta) given below refers 
to the center of the elliptical orbit of the eclipsing pair around the common center of mass 
of the triple, as opposed to the barycenter (see Irwin 1952). For consistency, corrections for 
light travel time in the inner binary have also been applied to the dates of the radial velocity 
observations at each iteration, which has only a minimal effect on the spectroscopic elements. 
These corrections range from —0.0145 days to ±0.0043 days (a change of 27 minutes over 
the 6.7-yr timespan of the velocities), and are listed in Table 1. 

Standard non-linear least-squares techniques were used to solve simultaneously for the 19 
unknowns that minimize the overall of the observations (Press et al. 1992). Weights were 
assigned to the measurements according to their individual errors. Most of the older times of 
minimum have no published errors, so those were determined by iterations in order to achieve 
a reduced near unity for each type of observation (visual, photographic, photoelectric). 
The adopted errors are 0.021 days for photographic data, 0.011 days for the visual timings, 
and 0.0015 days for photoelectric/CCD minima. The published errors of the more recent 
CCD timings were found to require adjustments to achieve reduced x^ values of unity. The 
scale factors used are 0.86 and 3.62 for the primary and secondary minima, respectively. 
The larger factor for the secondary is consistent with the smaller depth of that eclipse. 
The published uncertainties of the HIPPARCOS observations were found to require a minor 
adjustment by a factor of 0.85, and the optimal errors for the radial velocity measurements 
were determined to be 1.51 kms“^, 5.40 kms“^, and 1.34 kms“^ for stars Aa, Ab, and B, 
respectively. A total of 332 observations were used: 75 radial velocities for each star, 43 
times of primary eclipse, 6 times of secondary eclipse, 56 HIPPARCOS observations, and 2 
parallax measurements (see above). 

The results of this global £t are given in column 3 of Table 10, including derived quan¬ 
tities such as the masses of the three stars. The uncertainties listed for all derived quantities 
take full account of correlations between different elements (off-diagonal elements of the co- 
variance matrix). Initial solutions allowing for eccentricity in the inner orbit yielded a value 
not significantly different from zero, so a circular orbit was adopted. The outer orbit is 
moderately eccentric (e = 0.469), and its period is determined to be 15.8 yr, with an error 
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of only about 1.3%. 

The radial velocity observations and orbital solution for the inner eclipsing pair are 
shown graphically in Figure 9, with the motion in the outer orbit removed. In Figure 10a 
we show the velocities of the third star as a function of time, along with the measured 
velocities of the center of mass of the inner binary and the computed curves for the outer 
orbit. The 0 — C timing residuals (primary and secondary) from the best-£t linear ephemeris 
(Table 10) are displayed as a function of time in Figure 10b, for the more recent visual and 
photoelectric/CCD measurements. The solid line represents the computed third body effect. 
The semi-amplitude of this curve, as derived from our global fit, is Kq-c = 0.01458±0.00048 
days. Residuals from the curve for the individual timings are listed in Table 2. The time axis 
in both panels of Figure 10 is the same to show the complementarity of the data. Additionally, 
the time location of the HIPPARCOS observations is shown in the lower panel, and happens 
to be centered around periastron passage in the outer orbit. 

The astrometric motion of the photocenter of VlOOlCyg on the plane of the sky is 
illustrated in Figure 11, where the axes are parallel to the Right Ascension and Declination 
directions. The top panel shows the total motion resulting from the combined effects of 
annual parallax, proper motion, and orbital motion. The dominant contribution is from the 
proper motion (nearly 40 mas yr“^), which is indicated with an arrow. Parallax and orbital 
motion are smaller and comparable effects. In Figure 11b we have subtracted the proper mo¬ 
tion and parallactic contributions, leaving only the orbital motion with a 15.8-yr period and 
a semimajor axis of 9.6 mas. The direction of motion (retrograde) is indicated by the arrow. 
The individual HIPPARCOS observations are depicted schematically in both panels of this 
figure, but are seen more clearly in Figure 11b. Because they are one-dimensional in nature, 
the exact location of each measurement on the plane of the sky cannot be shown graphically. 
The Filed circles represent the predicted location on the computed orbit. The dotted lines 
connecting to each Filed circle indicate the scanning direction of the HIPPARCOS satellite 
for each measurement, and show which side of the orbit the residual is on. The short line 
segments at the end of and perpendicular to the dotted lines indicate the direction along 
which the actual observation lies, although the precise location is undetermined. Occasion¬ 
ally more than one measurement was taken along the same scanning direction, in which case 
two or more short line segments appear on the same dotted lines. Figure 11b shows that the 
observations are all on the long side of the orbit, and due to its orientation, they are roughly 
parallel to the direction of the proper motion. Since the orbit was not known at the time 
the HIPPARCOS catalog was released, we expect the orbital motion to have been absorbed 
into the published proper motion to some extent. Although the eFect is not large, we do see 
a hint of this in that both components of the proper motion as published are larger than the 
values we derive from our global solution, by about 2.5 mas yr“^ and 4.0 mas yr“^ for /i* 
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and ^jlsi respectively. Statistically the effect is somewhat more signihcant in Declination. 


4. Absolute dimensions 

The combination of the photometric resnlts in §2.3 and the spectroscopic resnlts in §3 
yields the absolnte masses and radii of the components of the eclipsing binary, which we list 
in Table 11 along with other properties as well as inferred parameters for the third star (see 
§5.3). For the photometric resnlts we have adopted the weighted average of the Gettysbnrg 
and URSA elements, which are rAa = 0.1668±0.0017, k = 0.598±0.004, and ix = 88?2±0?4. 
The precision of the absolnte masses of the binary components is 1.2% for star Aa and 0.7% 
for Ab, while that of star B is only 4%. The absolnte radii are determined to 1% and 2% 
for the primary and secondary, respectively. The measnred projected rotational velocities of 
both stars are in excellent agreement with the valnes compnted for synchronous rotation, as 
expected for such a short period. 

The parallax determination is much improved compared to the original HIPPARCOS 
result, the error being reduced from 17% to about 3% driven mostly by the individual 
estimates for stars Aa and Ab based on their physical properties, and to some extent also 
by the astrometry. The increased accuracy has an impact on the comparison with theory 
described below. The corresponding distance to the system is 166.9 ± 5.6 pc. The space 
motion of V1061Cyg in the Galactic frame is U = —20 kms“^, V = —2 kms“^, and 
W = +16 kms“^ (relative to the Local Standard of Rest), quite typical of Population I 
stars. The inclination angles of the inner and outer orbits (^a and ^ab) are known from 
our spectroscopic and light curve analyses, but the relative inclination 0 between the two 
orbits cannot be determined because the position angle of the node for the eclipsing pair 
is unknown^. However, a lower limit of 0niin = 20?2 ± 4?9 can be placed, which appears to 
exclude coplanarity. The semimajor axis of the photocenter of the binary in the wide orbit is 
only 9.55 mas, but the relative semimajor axis is much larger (55.3 mas). Gonsidering that 
the orbit is eccentric and therefore that the angular separation can be as large as 80 mas, the 
third star is potentially resolvable with current techniques (e.g., speckle interferometry on 4- 
m class telescopes, or adaptive optics on 8-10-m class telescopes) despite being 2 magnitudes 
fainter than the binary in the visible. Maximum separation should occur around the year 
2014. 


^The relative inclination angle of the two orbits is given by cos <j) = cos za cos iab + sin ia sin zab cos(Pa — 
Pab) (e.g., Fekel 1981). Since Pa is unknown, we can only set limits to cos(Pa — Pab) between —1 and +1, 
which leads to za — Gb < <5^' < G + Gb- 
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The central surface brightness parameter from the light curve solutions constrains the 
difference between the effective temperatures of the stars in the binary. To estimate the 
mean temperature of the system we made use of absolute UBV photometry outside of eclipse 
available from the work of Lacy (1992), the Bt and Vt magnitudes from the Tycho-2 catalog 
(Hpg et ah 2000), and the JHKg magnitudes from the 2MASS catalog. Color indices from 
various combinations of these magnitudes along with calibrations by Alonso, Arribas & 
Mart Inez-Roger (1996) and Ramirez & Melendez (2005) lead to a mean temperature close 
to 6000 K, similar to our estimate in §2.3. Folding in the spectroscopic estimates we obtain 
average values for the primary and secondary of 6180 ±100 K and 5300 ±150 K, respectively, 
which we adopt here. These correspond to spectral types of approximately F9 and G8, while 
the third star is roughly G4 (Cox 2000). 


5. Discussion 

The high precision of the mass and radius determinations for stars Aa and Ab (errors 
of 1-2%), along with the subsolar secondary mass and a mass ratio that is quite different 
from unity (g = 0.7266), make V1061Cyg a particularly valuable system to test models 
of single-star evolution. The four measured properties available for testing are the mass, 
radius, effective temperature, and luminosity (My). The latter two are of course related 
through the Stefan-Boltzmann law, although they were determined independently here since 
we know the distance. Because the temperatures were derived somewhat more indirectly and 
have greater uncertainties than My, we have preferred to rely here on the absolute visual 
magnitudes, inferred from the apparent brightness of the system (ignoring extinction; see 
§2.3), the fractional luminosities, and the parallax. 

No accurate estimate of the metallicity is available. The Population I kinematics (§4) 
provide only circumstancial evidence that the composition is perhaps near solar. In the 
following we have made that assumption initially, although non-solar compositions have also 
been explored by leaving the metallicity as a free parameter. In principle this additional 
freedom may appear to weaken the test; however, in this particular case it has no effect on 
our main conclusion, as seen below. 


5.1. Comparison with evolutionary models 

In Figure 12 we display the radius and absolute visual magnitude of V1061Gyg as a 
function of mass, along with theoretical isochrones from the Yonsei-Yale series of evolutionary 



models (Yi et al. 2001; Demarque et al. 2004) for solar metallicity (Zq = 0.01812). As 
do most models, these treat convection in the mixing-length approximation, and adopt a 
mixing-length parameter that best hts the properties of the Sun (ctml = 1.7432). The 
color transformations and bolometric corrections adopted are from the tables by Lejenne, 
Cuisinier & Buser (1998). The hgure shows that, within the errors, a 3.4-Gyr isochrone for 
this metallicity provides an excellent match to all observations except for the radius of the 
secondary, which is ~10% larger than predicted. This very signihcant difference is 5 times 
the nncertainty. Since the luminosity of the secondary is well reprodnced by the models, it 
follows that the effective temperature must be cooler than predicted. Indeed, our estimate 
from §4 is = 5300 ± 150 K, formally cooler than the predicted valne of 5460 K bnt only 
at the la level. A more precise estimate of the secondary temperature may be derived from 
our My and the radius of the star (with bolometric corrections taken also from Lejeune, 
Cnisinier & Buser 1998, for consistency). This gives = 5210 ±80 K, which is lower than 
the models predict at a more signihcant 3.1(T level. Thus the secondary does appear to be 
too cool by roughly 200 K. 

Extensive tests show that adjustments in the age or metallicity of the isochrones do not 
improve the agreement. Figure 13a illustrates this in the age/metallicity plane. The shaded 
areas represent all isochrones that are consistent with the measured mass, radius, and My of 
each star within the estimated uncertainties. The shaded areas do not overlap, meaning that 
no single isochrone can reproduce all properties of both stars simultaneously. The region 
spanned by models that are consistent with the primary properties ([Fe/H] within ±0.14 
dex of solar, ages of 3.06-3.78 Gyr) are shown in the mass/radius and mass/ My diagrams 
in Figures 13b and 13c. None of those models come close to matching the radius of the 
secondary star, although they do agree well with its absolute magnitude, as stated above. 

As a test we repeated the comparison using a different set of isochrones by Girardi et 
al. (2000), from the Padova group. In these models the adopted mixing-length parameter 
that best reprodnces the Snn (for Zq = 0.019) is ctml = 1.68, and a nnmber of other 
physical ingredients are somewhat different, as are the color transformations and bolometric 
corrections. Nevertheless, the hts to V1061Gyg are very similar, and again indicate a large 
discrepancy in the secondary radius R^h and a smaller difference in temperature. 

Both of the above sets of evolutionary models adopt a grey approximation for the 
outer boundary conditions to the internal stncture equations. More sophisticated non-grey 
atmospheres have been used in the models by Siess, Forestini & Dongados (1997) and Baraffe 
et al. (1998), which has been shown to be important for low mass stars (see, e.g., Ghabrier 
& Baraffe 1997). Additionally, these models incorporate an improved equation of state and 
other rehnements, although all this should only have an impact for stars considerably smaller 
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than the secondary of VlOhlCyg (which is near the calibration point at 1 M©). Indeed, 
the comparison between VlOblCyg and the Baraffe et al. (1998) isochrones {Zq = 0.020, 
c^ML = 1-9) is not very different from the two previous hts, and yields a similar age as before 
based on the agreement with all observed properties except i?Ab- 

The measured radii depend mainly on the results from the light curve fits and to a 
lesser degree on the spectroscopy. As described earlier, tests carried out during the analysis 
of both types of data make it very unlikely that either set of results is biased enough to 
explain the discrepancy. Since the properties of the more massive primary star appear to be 
well reproduced by the models (as expected from previous experience; see Andersen 2003), 
suspicion falls on our theoretical understanding of stars under a solar mass. Similar evidence 
has been presented over the last several years (e.g.. Popper 1997; Clausen et al. 1999a; Torres 
& Ribas 2002; Ribas 2003; Dawson & De Robertis 2004; Lopez-Morales & Ribas 2005), but 
in fact indications go as far back as the work by Hoxie (1973), Lacy (1977), and others. All 
of these studies have shown that theoretical calculations for stars less massive than the Sun 
tend to underestimate the radius by as much as 10% to 20%. Temperature differences also 
appear to be present, as in V1061 Cyg. 

The models by Baraffe et al. (1998) we used above are also available for lower values 
of the mixing length parameter in the mass range of interest. In Figure 14 we compare 
the observations for VlOOlCyg against solar metallicity models for auL = 1-0 and 1.5, in 
addition to q;ml = 1-9, adjusting the age to satisfy the constraints on the primary star 
within the errors. A 2.6-Gyr isochrone with ctml = 1-0 provides a signihcantly better £t to 
the secondary radius than the one with ctml = 1-9 used earlier, although the model prediction 
still falls short of the observed RAh value by about 1.8(T at the measured mass. The reduced 
convection implied by the lower mixing length parameter in these models compared to the 
previous ones has the effect of increasing the theoretical radii and slightly lowering the 
effective temperatures. The predicted temperature for ctml = 1-0 and M = MAh is 5110 K, 
somewhat closer than before to the empirical result of 5210 K (see above). 

The considerably smaller value of ctml apparently required by a star so close to the solar 
mass as VlOOlCyg Ab is perhaps somewhat surprising, although other studies have found 
similar evidence in some stars and have even argued for a mass dependence, with smaller 
values of ctml for later-type stars (e.g., Lastennet et al. 2003). However, the observational 
evidence for this is often contradictory (see Eggenberger et al. 2004). Theoretical studies, on 
the other hand, have tended to predict the opposite dependence with mass (Ludwig, Freytag 
& Steffen 1999; Trampedach et al. 1999). 



5.2. Comparison with other stars: the activity-radius connection 


Further clues on these disagreements may be gained from looking at other eclipsing 
binary systems with well determined properties. To avoid any possible dependence on mass, 
we focus here on binaries having at least one star in the same mass range as the secondary 
of VlOblCyg. Three such systems are available with reliable determinations of the mass, 
radius, and temperature good to better than 3%: RW Lac (Lacy, et ah 2005), in which 
the primary is within 0.5% of the mass of VlOOlCyg Ab; FL Lyr, with a secondary only 
3% more massive; and HS Aur, with a primary that is 3.5% less massive. The latter two 
are from work by Popper et al. (1986) and Andersen (1991), and FL Lyr has the added 
advantage that the components have quite dissimilar masses (as in the case of V1061Cyg), 
which provides increased leverage for the test. The main properties of these systems are 
collected in Table 12, along with V1061Cyg and the Sun for reference. Two additional 
systems (CG Cyg and RT And; Popper 1994) have one component in the same mass range, 
but the temperature determinations are considerably more uncertain and are thus less useful 
for our purposes. 

We compared the measurements for RW Lac, FL Lyr, and HS Aur to the Yonsei-Yale 
models in the same way we did for V1061 Cyg, adjusting the age and metallicity of the 
isochrones to obtain the best possible £t (see Figure 15). FL Lyr is seen to present the same 
problem as V1061Cyg: the secondary appears too large compared to theory, by about the 
same amount as we saw before (~10%), while the properties of primary are well fit. On the 
other hand, the primary of RW Lac is virtually identical in mass to V1061Cyg Ab yet it 
shows no indication of a radius discrepancy. The primary in HS Aur also seems to be well 
reproduced by the models (as well as the secondary), within the errors. 

A pattern that may explain why stars of very similar mass sometimes appear too large, 
while other times they conform well to theory, is seen in the activity level they present. 
Both V1061Cyg Ab and FL Lyr B, which show the radius discrepancy, are in relatively 
tight systems with orbital periods of only 2.35 days and 2.18 days, respectively. The stars 
in these binaries are rapid rotators (consistent with synchronous rotation maintained by 
tidal forces): V1061Cyg Ab has a measured usini = 20 ± 3 kms“^ (§2.1) and FL Lyr B 
has usini = 25 ± 2 kms“^ (Popper et al. 1986). Both binaries are strong ROSAT X-ray 
sources (Voges et al. 1999) and their X-ray luminosities (Table 12) are in line with those of 
active single and binary stars with similar usini values (e.g., Cutispoto et al. 2003, Fig. 5). 
Although to our knowledge no observations exist to verify whether the Ca II H and K lines 
are in emission, both systems show signs of intrinsic variability in the light curve suggesting 
the presence of spots (see Popper et al. 1986 for FL Lyr, and §2.3 for V1061Cyg). They 
are thus active binaries. RW Lac and HS Aur, on the other hand, do not show a radius 
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discrepancy and happen to have mnch longer orbital periods of 10.37 days and 9.82 days, 
respectively, and their components are slow rotators. No Ca 11 emission is seen in HS Anr 
(Popper 1976). Neither binary was detected as an X-ray source by ROSAT (Voges et al. 
1999), and by all accounts they appear inactive. A direct relation is thus seen for stars of the 
same mass between the activity level and the increased stellar size compared to predictions 
from standard models (i.e., those adopting a mixing length parameter matching the Sun): 
active stars are larger, and inactive ones appear normal. The evidence for V1061Cyg and 
other later-type stars also indicates that active stars are cooler. 

This seems to have been mentioned only as a possibility (among several others) in 
previous studies reporting radius discrepancies for stars under a solar mass, at least those 
originating from the eclipsing binary community. But as a matter of fact, the connection 
between activity and the global properties of low mass stars has been studied in some detail 
previously in a slightly different context, although without direct and precise knowledge of the 
masses and radii of the stars involved. Mullan & MacDonald (2001) investigated the effects 
of magnetic fields on the sizes and effective temperatures of active versus inactive M dwarfs, 
and found empirical evidence that a higher activity level leads to larger radii and cooler tem¬ 
peratures. Their sample consisted of single stars with effective temperature determinations 
from infrared spectroscopy and bolometric luminosities from multi-band photometry, from 
which stellar radii were inferred indirectly with typical uncertainties of 10-15%. Although 
their work focussed mainly on the consequences for the internal structure of fully convective 
stars, their initial attempts at modeling magnetic helds were successful in describing these 
effects to hrst order. The present study on V1061Cyg shows the activity-radius connection 
clearly for stars with accurately determined dynamical masses, radii, and temperatures good 
to 1-3%. More importantly, the effect is seen for objects that are only 7% less massive than 
the Sun. 

A theoretical understanding of this connection is also not new, and provides some insight 
into the better agreement between the observations for V1061 Cyg and the 1ow-q;ml models 
in Figure 14. Strong magnetic helds commonly associated with chromospheric activity have 
been shown to inhibit the efficiency of convective heat transport (e.g., Bray & Loughhead 
1964; Gough & Tayler 1966; Stein, Brandenburg & Nordlund 1992, and references therein), 
and as a result the size of the star must grow larger to radiate away the same amount of 
energy. The decreased convection effectively leads to a lower value of the mixing length 
parameter (see Tayler 1987), which explains the better £t to the radius of the secondary 
using the Baraffe et al. (1998) models with uml = 1-0 as opposed to those with the solar 
value of ctml = 1-9- Similar improvements in the fit to other low-mass stars using a reduced 
mixing length parameter were also reported by Clausen et al. (1999a), and in fact much 
earlier by Gabriel (1969) and Cox, Shaviv & Hodson (1981). 
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The enhanced activity in VlOOlCyg Ab is driven by rapid (synchronous) rotation in¬ 
duced by tidal forces. The primary star is rotating even more rapidly (n sini = 36±2 kms“^), 
yet it shows no obvious indication of a signihcantly larger radius compared to standard mod¬ 
els. This is most likely because it is a more massive star (spectral type F9, M = 1.282 Mq) 
and therefore its convective envelope (where magnetic activity takes place) is signihcantly 
reduced. To illustrate this in a more quantitative way, Figure 16 shows the mass of the 
convective envelope for main-sequence stars as a function of stellar mass from the models 
by Siess, Forestini & Dougados (1997) (see also Siess, Dufour & Forestini 2000), for solar 
composition and the age of 3.4 Gyr we infer for the system. In the lower panel the envelope 
mass is expressed as a percentage of the stellar mass. The inset shows an enlargement of 
the region relevant to V1061 Cyg. According to these models the convective envelope of the 
secondary represents about 4.7% of its total mass, whereas that of the primary is only 0.2% 
(~20 times smaller). A similarly reduced convective envelope for the primary of FL Lyr 
(M = 1.221 Mq) explains why that star also does not show evidence of an enlarged radius, 
while the secondary (which is similar in mass to V1061 Cyg Ab) does. 

The extent to which a star is enlarged by this effect may be expected to depend upon 
the strength of the activity. We note, for instance, that although models with oml = 1-0 
seem to match the radius of V1061 Cyg Ab better, the same models still do not reproduce 
the properties of apparently more active lower-mass stars such as YY Gem (Torres & Ribas 
2002), CU Cue (Ribas 2003), GU Boo (Lopez-Morales & Ribas 2005), and others. It is clear 
that further examples of stars (both active and inactive) in the suitable mass range with 
accurately determined parameters are needed to explore this. Systems such as V1061Cyg 
and FL Lyr, with components of appreciably different mass and primaries that are larger 
than ~1.2 Mq, are particularly useful because they allow the activity-radius effect to be 
separated out, since the primary should not be signihcantly affected. 


5.3. The tertiary star 

The V1061Cyg system is especially valuable for containing a third star of mass in¬ 
distinguishable from that of the active star Ab, but with slow rotation. It is presumably 
an inactive analog of the secondary, with the same age and chemical composition, and is 
thus ideal for testing the effect discussed above. We do not have a direct measurement of 
its size, but the radius can be inferred from the effective temperature and luminosity, and 
is Rb = 0.870 ± 0.087 R©. While this is formally smaller than the secondary radius of 
Ras = 0.974 ±0.020 R©, as we expected, it is not precise enough for a dehnitive test. Other 
inferred properties of the tertiary are listed in Table 11. The star is represented graphically 
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in Figure 15, where it is seen to be in good agreement with the isochrone for the primary, 
within the admittedly large errors. The effective temperature we determine for the tertiary 
(5670 ± 150 K) is somewhat hotter than that of the secondary even though the masses are 
similar, which again is consistent with our conclusion that the activity in the secondary has 
made that star cooler. 

The formal uncertainty in is dominated by errors in the light fraction £b and effective 
temperature Tg^, in that order of importance. Improvements in the former could be made, 
for example, by direct detection of the third star and measurement of its brightness through 
speckle interferometry or adaptive optics imaging (§4). A better temperature estimate would 
require spectroscopy with higher S/N, although the triple-lined nature of the spectrum and 
the faintness of the star would still pose a challenge. Additional radial velocity measurements 
over the coming years covering the next nodal passage in the outer orbit (in 2009.1) would 
also be important to refine the value of the tertiary mass. 


6. Concluding remarks 

The results of our spectroscopic, photometric, and astrometric analyses of the V1061 Cyg 
system have taken us in a rather different direction than we anticipated when we began this 
study. The possible status of the object as an example of the rare class of cool Algols is now 
clearly ruled out. Instead, we have shown here that it is a hierarchical triple system with 
an outer period of 15.8 yr, in which the eclipsing inner pair is well detached, is composed of 
main-sequence stars, has a mass ratio quite different from unity, and has a secondary that 
is slightly below a solar mass. The absolute masses and radii for the binary components are 
determined with a relative precision of 2% or better, and the mass of the third star is good 
to 4%. 

While the primary star is well fit in mass, radius, temperature, and luminosity by 
standard stellar evolution models with a metallicity near solar and a mixing length parameter 
set by the calibration to the Sun, the secondary appears ~10% too large. This discrepancy 
is 5 times the size of the observational errors, and quite surprising for a star that differs by 
only 7% in mass from the Sun. There are also indications that it is cooler than predicted by 
some 200 K. V1061Cyg is yet another example highlighting our incomplete understanding 
of the structure and evolution of stars in the lower main sequence. Similar differences in size 
and temperature have been noticed previously for lower-mass stars with accurately measured 
properties, but the source of the problem has remained unclear. By comparing V1061Cyg 
Ab to several other objects of nearly the same mass we have identified the activity level 
as a key factor distinguishing cases that show the radius discrepancy from those that do 
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not. This link between activity and increased radius has been mentioned in the literature 
before, but is shown here for the first time for stars with accurately known masses, radii, 
and temperatures. 

It is often stated that the structure and evolution of stars (particularly those close to a 
solar mass) are completely determined once the chemical composition and mass are specified. 
It is quite clear now that for stars of order 1 Mq or less an additional parameter must be 
taken into account, which has to do with the level of chromospheric activity. Whether this 
parameter is directly the rotational speed (e.g., nsini) or period, the magnetic field strength, 
the Rossby number (e.g., Noyes et ah 1984; Basri 1987), or some other more complicated 
activity indicator remains to be determined. To first order it appears that the effective mixing 
length may be a useful proxy, but with one exception current stellar evolution models are 
not publicly available for more than one value of oml, so testing this is somewhat difficult 
in practice. If the predictions from the models are to reach an accuracy matching current 
observations of low-mass stars (~l-2% relative errors in the masses and radii), this effect 
can no longer be ignored. Further progress will require more examples of binary systems 
with well determined physical parameters and different levels of activity in the relevant mass 
regime in order to help calibrate any such parameter. 

Testing models of single-star evolution by means of eclipsing binaries, as astronomers 
have done for decades, might perhaps be seen as part of the problem since the enhanced 
activity displayed by many of these systems is a direct result of tidal synchronization that 
occurs only in close binaries. Although wider eclipsing pairs with inactive components under 
1 Mq certainly do exist, they are less common and in some respects more difficult to study. 
Stars near the bottom of the main sequence, M dwarfs in particular, are found to be active 
more often than not. 

The effects of magnetic fields on the evolution of stars have already begun to be ex¬ 
plored by theorists (e.g., D’Antona, Ventura & Mazzitelli 2000; Mullan & MacDonald 2001; 
Maeder & Meynet 2003, although the latter authors focus on more massive stars), and ini¬ 
tial comparisons with observations are encouraging. One area where this is likely to have 
a significant impact is the study of T Tauri stars, which are typically very active. Even in 
substellar objects such as brown dwarfs activity appears to be quite common, and might be 
expected to have similar consequences on their structure. Some evidence for this has already 
been reported (Mohanty, Jayawardhana & Basri 2004). 
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Table 1. Radial velocity measurements for VlOGlCyg. 


HJD 

(2,400,000+) 

Year 

RVAa 

(kms“^) 

RVAb 
(km s“^) 

RVb 
( km s“^) 

(O-C)Aa 

(kms“^) 

(O-C)Ab 

(kms“^) 

(O-C)b 

(kms~^) 

Inner 

Phase^ 

Outer 

Phase^ 

AT*" 

(days) 

47101.6142 

1987.8347 

+66.13 

-106.02 

+2.24 

-0.13 

+7.64 

-1.13 

0.8346 

0.7504 

-0.0045 

51039.8277 

1998.6169 

-35.07 

+30.15 

-7.37 

-0.30 

-5.64 

-0.27 

0.0549 

0.4313 

-0.0143 

51057.7385 

1998.6660 

+74.19 

-119.55 

-8.69 

-1.80 

-2.80 

-1.69 

0.6874 

0.4343 

-0.0143 

51075.7492 

1998.7153 

-71.66 

+84.41 

-5.67 

+0.32 

-2.39 

+1.22 

0.3624 

0.4375 

-0.0144 

51087.6477 

1998.7478 

-40.88 

+42.78 

-5.91 

+0.30 

-1.56 

+0.91 

0.4328 

0.4395 

-0.0144 

51089.6402 

1998.7533 

-90.48 

+115.25 

-6.76 

+0.79 

+1.99 

+0.05 

0.2819 

0.4399 

-0.0144 

51092.6176 

1998.7614 

+22.87 

-41.20 

-5.72 

+0.59 

+1.83 

+1.07 

0.5507 

0.4404 

-0.0144 

51114.6251 

1998.8217 

+31.92 

-58.71 

-4.83 

-0.77 

-1.23 

+1.84 

0.9289 

0.4442 

-0.0144 

51122.6356 

1998.8436 

-77.41 

+105.84 

-5.94 

+1.28 

+10.08 

+0.68 

0.3425 

0.4456 

-0.0144 

51146.6110 

1998.9093 

+25.53 

-40.99 

-6.31 

-1.10 

+8.35 

+0.17 

0.5593 

0.4497 

-0.0145 

51163.5148 

1998.9556 

+83.04 

-121.98 

-5.20 

+0.87 

+3.89 

+1.18 

0.7627 

0.4526 

-0.0145 

51356.6644 

1999.4844 

-42.89 

+46.35 

-4.87 

+0.90 

-0.04 

+0.40 

0.0711 

0.4860 

-0.0145 

51423.5273 

1999.6674 

+28.65 

-48.61 

-5.95 

+0.28 

+4.70 

-1.06 

0.5640 

0.4976 

-0.0145 

51464.6818 

1999.7801 

-56.90 

+74.95 

-5.01 

+1.50 

+9.06 

-0.35 

0.1015 

0.5047 

-0.0144 

51473.6526 

1999.8047 

+32.28 

-64.30 

-5.45 

-1.79 

-2.87 

-0.84 

0.9244 

0.5063 

-0.0144 

51477.5968 

1999.8155 

+48.61 

-77.00 

-3.19 

+0.85 

+3.29 

+1.40 

0.6051 

0.5069 

-0.0144 

51510.5931 

1999.9058 

+69.96 

-115.42 

-4.44 

+0.25 

-4.73 

-0.04 

0.6661 

0.5126 

-0.0143 

51537.4847 

1999.9794 

-68.89 

+85.11 

-3.12 

-0.24 

+5.52 

+1.13 

0.1257 

0.5173 

-0.0143 

51558.4542 

2000.0368 

-38.41 

+40.50 

-3.26 

+1.08 

+1.16 

+0.87 

0.0616 

0.5209 

-0.0143 

51641.8900 

2000.2653 

+51.05 

-85.08 

-2.42 

-1.32 

+2.47 

+1.24 

0.6169 

0.5353 

-0.0140 

51713.6479 

2000.4617 

-89.95 

+99.16 

-3.43 

-0.48 

-8.10 

-0.17 

0.1958 

0.5477 

-0.0138 

51715.6836 

2000.4673 

-40.85 

+37.33 

-2.04 

-0.13 

-2.83 

+1.21 

0.0633 

0.5481 

-0.0138 

51777.8307 

2000.6375 

+14.28 

-33.17 

-3.10 

-4.29 

+8.63 

-0.20 

0.5467 

0.5588 

-0.0136 

51797.8286 

2000.6922 

-43.62 

+35.79 

-2.25 

-0.04 

-7.84 

+0.54 

0.0686 

0.5623 

-0.0135 

51842.6036 

2000.8148 

-77.47 

+94.87 

-1.64 

+0.26 

+4.49 

+0.90 

0.1490 

0.5700 

-0.0133 

51857.6308 

2000.8559 

+21.79 

-42.67 

-1.56 

+0.23 

+3.68 

+0.89 

0.5527 

0.5726 

-0.0133 

51861.5475 

2000.8667 

-93.89 

+112.74 

-3.22 

-0.41 

+0.79 

-0.79 

0.2218 

0.5733 

-0.0132 

51888.5267 

2000.9405 

+79.34 

-120.56 

-3.76 

+0.30 

+5.08 

-1.48 

0.7187 

0.5780 

-0.0131 

51902.4960 

2000.9788 

+70.32 

-112.81 

-2.21 

+0.06 

+0.83 

-0.01 

0.6716 

0.5804 

-0.0130 

51921.4546 

2001.0307 

+80.78 

-129.06 

-2.76 

+0.13 

-1.02 

-0.67 

0.7507 

0.5837 

-0.0130 

52071.7515 

2001.4422 

+78.34 

-115.61 

-1.55 

+2.05 

+7.26 

-0.30 

0.7983 

0.6097 

-0.0121 

52099.9064 

2001.5193 

+76.42 

-117.95 

-1.80 

-0.13 

+5.44 

-0.71 

0.7963 

0.6145 

-0.0119 

52105.5866 

2001.5348 

-92.25 

+104.57 

-0.43 

+1.28 

-6.09 

+0.63 

0.2168 

0.6155 

-0.0119 

52121.6122 

2001.5787 

-31.49 

+36.17 

+1.78 

+1.19 

+9.35 

+2.75 

0.0460 

0.6183 

-0.0118 

52158.5190 

2001.6797 

+77.10 

-129.27 

-0.11 

-2.04 

-1.99 

+0.65 

0.7735 

0.6247 

-0.0116 

52182.6433 

2001.7458 

-36.58 

+34.81 

-0.92 

+0.38 

+2.44 

-0.29 

0.0539 

0.6288 

-0.0114 

52194.6505 

2001.7786 

-83.42 

+109.65 

-0.30 

+1.53 

+11.30 

+0.26 

0.1707 

0.6309 

-0.0113 

52220.6532 

2001.8498 

-97.12 

+112.70 

-0.65 

-1.43 

-0.30 

-0.24 

0.2515 

0.6354 

-0.0111 

52226.5786 

2001.8661 

+78.36 

-127.09 

-0.74 

-0.36 

-0.01 

-0.36 

0.7766 

0.6364 

-0.0111 

52234.5686 

2001.8879 

-87.85 

+98.85 

-1.62 

-0.15 

-3.07 

-1.29 

0.1814 

0.6378 

-0.0110 

52240.5913 

2001.9044 

+77.39 

-128.49 

+1.00 

-2.51 

+0.29 

+1.30 

0.7479 

0.6388 

-0.0110 

52247.5736 

2001.9235 

+79.01 

-123.51 

-0.23 

+0.35 

+3.62 

+0.03 

0.7234 

0.6401 

-0.0109 

52266.4777 

2001.9753 

+78.62 

-130.69 

-1.33 

+0.25 

-3.87 

-1.18 

0.7792 

0.6433 

-0.0108 

52274.4698 

2001.9972 

-91.28 

+103.79 

-1.42 

-2.68 

+0.86 

-1.31 

0.1850 

0.6447 

-0.0107 

52287.4668 

2002.0328 

+78.79 

-123.26 

-1.55 

+0.20 

+3.98 

-1.52 

0.7236 

0.6470 

-0.0106 
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HJD 

(2,400,000+) 

Year 

RVAa 

(kms“^) 

RVAb 
(km s“^) 

RVb 
( km s“^) 

(O-C)Aa 

(kms~^) 

(O-C)Ab 

(kms~^) 

(O-C)b 

(kms~^) 

Inner 

Phase^ 

Outer 

Phase^ 

AT*" 

(days) 

52416.7769 

2002.3868 

+67.27 

-115.52 

+0.12 

-1.87 

-0.55 

-0.58 

0.8280 

0.6693 

-0.0095 

52446.6840 

2002.4687 

+29.45 

-61.17 

-0.44 

-0.89 

+0.56 

-1.31 

0.5727 

0.6745 

-0.0092 

52458.6699 

2002.5015 

+72.05 

-117.36 

+3.17 

+0.92 

+0.58 

+2.23 

0.6804 

0.6766 

-0.0091 

52537.6959 

2002.7179 

-78.33 

+90.48 

+2.81 

-0.91 

+4.42 

+1.42 

0.3568 

0.6902 

-0.0084 

52538.7321 

2002.7207 

+76.54 

-122.17 

-0.43 

+1.36 

+1.80 

-1.82 

0.7983 

0.6904 

-0.0084 

52566.6749 

2002.7972 

+74.98 

-123.61 

+2.60 

-0.81 

+1.36 

+1.05 

0.7060 

0.6952 

-0.0081 

52613.5082 

2002.9254 

+66.53 

-112.28 

+1.19 

+0.13 

+0.03 

-0.63 

0.6636 

0.7033 

-0.0076 

52626.4669 

2002.9609 

-89.56 

+100.34 

+1.01 

+0.07 

-2.02 

-0.88 

0.1859 

0.7056 

-0.0075 

52638.5178 

2002.9939 

-89.41 

+94.18 

+1.77 

-1.39 

-5.90 

-0.19 

0.3213 

0.7076 

-0.0074 

52651.4778 

2003.0294 

+63.56 

-112.08 

+1.00 

-0.46 

-2.84 

-1.03 

0.8441 

0.7099 

-0.0072 

52750.8421 

2003.3014 

-88.59 

+105.04 

+2.82 

+1.69 

+2.49 

+0.22 

0.1875 

0.7271 

-0.0061 

52788.7918 

2003.4053 

-77.02 

+92.33 

+6.46 

+0.03 

+8.21 

+3.65 

0.3596 

0.7336 

-0.0057 

52797.7304 

2003.4298 

-89.02 

+96.14 

+5.91 

-3.17 

-0.06 

+3.05 

0.1687 

0.7352 

-0.0056 

52891.7370 

2003.6872 

-101.07 

+113.83 

+2.46 

-4.54 

+3.47 

-0.93 

0.2290 

0.7514 

-0.0044 

52912.7271 

2003.7446 

-88.21 

+100.85 

+3.76 

-0.74 

+3.07 

+0.25 

0.1738 

0.7551 

-0.0041 

52943.5978 

2003.8292 

-86.20 

+104.52 

+3.30 

+0.56 

+7.89 

-0.38 

0.3292 

0.7604 

-0.0037 

52978.5239 

2003.9248 

-95.01 

+106.13 

+2.53 

+0.09 

-1.78 

-1.35 

0.2127 

0.7664 

-0.0033 

53006.4893 

2004.0013 

-74.36 

+79.73 

+4.33 

-0.60 

+1.33 

+0.30 

0.1300 

0.7713 

-0.0029 

53030.4508 

2004.0669 

-83.70 

+91.83 

+4.00 

-0.09 

+0.00 

-0.16 

0.3411 

0.7754 

-0.0026 

53127.8133 

2004.3335 

+67.91 

-114.52 

+4.48 

+1.38 

+0.82 

-0.20 

0.8316 

0.7922 

-0.0012 

53181.6702 

2004.4810 

+75.66 

-131.09 

+3.90 

-0.23 

-2.60 

-1.05 

0.7824 

0.8015 

-0.0004 

53187.6941 

2004.4975 

-81.15 

+92.63 

+6.52 

+0.20 

+4.74 

+1.54 

0.3495 

0.8026 

-0.0003 

53202.6395 

2004.5384 

+74.82 

-131.70 

+4.91 

-1.12 

-3.04 

-0.15 

0.7184 

0.8052 

-0.0001 

53275.6923 

2004.7384 

+62.31 

-103.54 

+4.43 

+1.39 

+4.79 

-0.98 

0.8495 

0.8178 

+0.0010 

53304.6418 

2004.8176 

-93.92 

+88.17 

+5.98 

-2.70 

-12.75 

+0.44 

0.1862 

0.8228 

+0.0015 

53323.6007 

2004.8695 

-98.17 

+111.57 

+5.96 

-0.36 

+1.66 

+0.34 

0.2654 

0.8261 

+0.0018 

53342.5272 

2004.9214 

-86.33 

+102.44 

+3.49 

+0.84 

+7.26 

-2.21 

0.3308 

0.8294 

+0.0021 

53357.4967 

2004.9623 

+72.33 

-127.89 

+5.58 

-2.28 

-0.35 

-0.18 

0.7100 

0.8319 

+0.0023 

53372.4574 

2005.0033 

-54.77 

+49.58 

+4.99 

+0.63 

-1.77 

-0.83 

0.0855 

0.8345 

+0.0025 

53480.8610 

2005.3001 

-96.41 

+112.01 

+8.87 

+0.37 

+4.08 

+2.68 

0.2812 

0.8533 

+0.0043 


^Computed from the ephemeris in the last column of Table 10. 
^Corrections for light travel time (see text). 








Table 2. Measured times of eclipse for V1061 Cyg. 


HJD 

(2,400,000+) 

a 

(days) 

Year 

Type*^ 

Eclipse'^ E 

0-C 

(days) 

Ref 

26355.245 

0.021 

1931.0342 

Pg 

1 

-10935.0 

+0.02093 

1 

26512.460 

0.021 

1931.4646 

Pg 

1 

-10868.0 

+0.00762 

1 

26925.476 

0.021 

1932.5954 

pg 

1 

-10692.0 

+0.00714 

1 

26958.310 

0.021 

1932.6853 

pg 

1 

-10678.0 

-0.01237 

1 

26958.335 

0.021 

1932.6854 

pg 

1 

-10678.0 

+0.01263 

1 

27183.594 

0.021 

1933.3021 

pg 

1 

-10582.0 

-0.00930 

1 

28692.483 

0.021 

1937.4332 

pg 

1 

-9939.0 

-0.02197 

1 

28753.487 

0.021 

1937.6002 

pg 

1 

-9913.0 

-0.03072 

1 

36788.465 

0.021 

1959.5988 

pg 

1 

-6489.0 

+0.02863 

1 

36868.261 

0.021 

1959.8173 

pg 

1 

-6455.0 

+0.03879 

1 

44166.347 

0.011 

1979.7983 

V 

1 

-3345.0 

+0.00712 

2 

45229.383 

0.011 

1982.7088 

V 

1 

-2892.0 

+0.00127 

2 

45229.385 

0.021 

1982.7088 

pg 

1 

-2892.0 

+0.00327 

2 

45933.377 

0.011 

1984.6362 

V 

1 

-2592.0 

-0.00071 

2 

46001.438 

0.011 

1984.8226 

V 

1 

-2563.0 

+0.00756 

2 

46651.4515 

0.0015 

1986.6022 

pe 

1 

-2286.0 

+0.00191 

3 

46705.420 

0.011 

1986.7500 

V 

1 

-2263.0 

-0.00215 

2 

47355.442 

0.011 

1988.5296 

V 

1 

-1986.0 

+0.00431 

2 

47362.476 

0.0015 

1988.5489 

pe 

1 

-1983.0 

-0.00095 

2 

48972.2716 

0.0007 

1992.9563 

pe 

1 

-1297.0 

-0.00127 

2 

48972.2735 

0.0005 

1992.9563 

pe 

1 

-1297.0 

+0.00063 

2 

49528.413 

0.011 

1994.4789 

V 

1 

-1060.0 

-0.02707 

2 

49535.4799 

0.0005 

1994.4982 

pe 

1 

-1057.0 

-0.00024 

4 

49535.4802 

0.0008 

1994.4982 

pe 

1 

-1057.0 

+0.00006 

4 

49941.4566 

0.0003 

1995.6097 

pe 

1 

-884.0 

-0.00046 

5 

49941.4572 

0.0003 

1995.6097 

pe 

1 

-884.0 

+0.00014 

5 

50286.4260 

0.011 

1996.5542 

V 

1 

-737.0 

+0.00719 

2 

51159.3771 

0.0022 

1998.9442 

pe 

1 

-365.0 

-0.00136 

2 

51159.3787 

0.0008 

1998.9442 

pe 

1 

-365.0 

+0.00024 

2 

51159.3789 

0.0009 

1998.9442 

pe 

1 

-365.0 

+0.00044 

6 

51159.3791 

0.0007 

1998.9442 

pe 

1 

-365.0 

+0.00064 

2 

52015.90554 

0.00010 

2001.2893 

pe 

1 

0.0 

+0.00013 

7 






32 


Table 2—Continued 


HJD 

(2,400,000+) 

cr 

(days) 

Year 

Type'^ 

Eclipse’’ 

E 

0-C 

(days) 

Ref 

52095.6909 

0.0002 

2001.5077 

pe 

1 

+34.0 

-0.00029 

8 

52102.73118 

0.00013 

2001.5270 

pe 

1 

+37.0 

+0.00008 

8 

52109.77080 

0.00012 

2001.5463 

pe 

1 

+40.0 

-0.00022 

8 

52149.6636 

0.0003 

2001.6555 

pe 

1 

+57.0 

-0.00029 

8 

52448.8603 

0.0011 

2002.4746 

pe 

2 

+184.5 

+0.00031 

9 

52482.8861 

0.0003 

2002.5678 

pe 

1 

+199.0 

-0.00006 

9 

52589.6558 

0.0018 

2002.8601 

pe 

2 

+244.5 

-0.00210 

9 

52602.56431 

0.00010 

2002.8955 

pe 

1 

+250.0 

-0.00006 

9 

52609.60438 

0.00010 

2002.9147 

pe 

1 

+253.0 

+0.00012 

9 

52786.7736 

0.0014 

2003.3998 

pe 

2 

+328.5 

-0.00113 

10 

52813.7610 

0.0002 

2003.4737 

pe 

1 

+340.0 

+0.00007 

10 

52834.8804 

0.0003 

2003.5315 

pe 

1 

+349.0 

-0.00017 

10 

52867.7336 

0.0003 

2003.6214 

pe 

1 

+363.0 

+0.00027 

10 

52887.6813 

0.0014 

2003.6761 

pe 

2 

+371.5 

+0.00165 

10 

52907.62597 

0.00012 

2003.7307 

pe 

1 

+380.0 

+0.00001 

10 

53145.8083 

0.0022 

2004.3828 

pe 

2 

+481.5 

+0.00014 

11 

53293.6434 

0.0011 

2004.7875 

pe 

2 

+544.5 

-0.00178 

11 


“^Technique: ph = photographic; v = visual; pe = photoelectric/CCD. 
'’I = primary; 2 = secondary. 


References. — 1. Strohmeier, Knigge & Ott (1962); 2. Kreiner, Kim & Nha 
(2000); 3. Diethelm (1986); 4. Agerer & Hiibscher (1995); 5. Agerer &: Hiibscher 
(1996); 6. Ogloza, Drozdz & Zola (2000); 7. Lacy (2001); 8. Lacy, Straughn & 
Denger (2002); 9. Lacy (2002); 10. Lacy (2003); 11. Lacy (2004). 


Note. — Timing uncertainties (a) have been determined or adjusted as described 
in the text. O—C residuals are from the global solution in §3. 








Table 3. 


Differential D-band photometry for V1061 Cyg from the URSA telescope. 


HJD 

(2,400,000+) 

AU 

(mag) 

52102.73110 

-0.101 

52428.91322 

-0.091 

52095.69134 

-0.091 

52076.91832 

-0.083 

52149.66397 

-0.102 


Note. — Table 3 
is published in its en¬ 
tirety in the electronic 
edition of the Astrophys- 
ical Journal. A portion 
is shown here for guid¬ 
ance regarding its form 
and content. 



Table 4. 


Differential 5-band photometry for V1061 Cyg from Gettysburg. 


HJD 

(2,400,000+) 

AB 

(mag) 

52548.59272 

0.896 

52541.55330 

0.907 

52548.59492 

0.892 

52548.59712 

0.885 

52541.55930 

0.900 


Note. — Table 4 is 
published in its entirety 
in the electronic edi¬ 
tion of the Astrophysi- 
cal Journal. A portion 
is shown here for guid¬ 
ance regarding its form 
and content. 



Table 5. 


Differential D-band photometry for V1061 Cyg from Gettysburg. 


HJD 

(2,400,000+) 

AD 

(mag) 

52541.55170 

0.592 

52548.59332 

0.575 

52541.55390 

0.590 

52548.59552 

0.577 

52548.59782 

0.575 


Note. — Table 5 is 
published in its entirety 
in the electronic edi¬ 
tion of the Astrophysi- 
cal Journal. A portion 
is shown here for guid¬ 
ance regarding its form 
and content. 



Table 6. 


Differential -R-band photometry for VI061 Cyg from Gettysburg. 


HJD 

(2,400,000+) 

AR 

(mag) 

52541.55210 

0.389 

52548.59372 

0.377 

52541.55430 

0.392 

52548.59592 

0.368 

52548.59822 

0.376 


Note. — Table 6 is 
published in its entirety 
in the electronic edi¬ 
tion of the Astrophysi- 
cal Journal. A portion 
is shown here for guid¬ 
ance regarding its form 
and content. 



Table 7. 


Differential /-band photometry for VI061 Cyg from Gettysburg. 


HJD 

(2,400,000+) 

AI 

(mag) 

52548.59172 

0.191 

52541.55240 

0.187 

52548.59402 

0.190 

52541.55460 

0.172 

52548.59622 

0.185 


Note. — Table 7 is 
published in its entirety 
in the electronic edi¬ 
tion of the Astrophysi- 
cal Journal. A portion 
is shown here for guid¬ 
ance regarding its form 
and content. 



Table 8. 


Comparison stars for V1061 Cyg. 


Star 

R.A. (J2000) 

Dec. (J2000) 

Rj' 

Vt 

URSA 

CSC 03600-00278 

21:07:34.082 

+52:06:51.60 

11.210 

10.605 

CSC 03600-00377 

21:07:22.127 

+52:01:37.07 

12.555 

11.824 

CSC 03600-00423 

21:06:55.420 

+52:02:49.14 

11.837 

11.456 

CSC 03600-00443 

21:06:40.440 

+52:08:02.96 

13.087 

11.451 

Gettysburg 

CSC 03600-00425 

21:06:49.626 

+52:00:17.08 

9.261 

9.026 

CSC 03600-00259 

21:07:40.472 

+51:52:56.94 

11.610 

10.318 


Note. — Coordinates and magnitudes and Vt are taken from 
the Tycho-2 Catalog (Hpg et al. 2000). 







Table 9. 


Light curve solutions for V1061 Cyg. 


Parameter^ 


Gettysburg 



URSA 

V 

B 

V 

R 

I 

Mean 

J = c/al/'^A a • • • • 

0.398 

0.472 

0.525 

0.590 


0.436 


0.005 

0.005 

0.006 

0.005 


0.002 

f'Aa. 

0.1635 

0.1625 

0.1635 

0.1648 

0.164 

0.1674 


0.0009 

0.0010 

0.0010 

0.0009 

0.002 

0.0009 

f'Ab. 

0.0959 

0.0970 

0.0971 

0.0994 

0.097 

0.1009 


0.0007 

0.0008 

0.0008 

0.0007 

0.002 

0.0006 

k = 7-Ab/»’Aa. 

0.587 

0.597 

0.594 

0.603 

0.595 

0.603 


0.003 

0.003 

0.003 

0.002 

0.007 

0.009 

lA (deg) . 

89.1 

88.5 

88.5 

88.3 

88.6 

87.9 


0.5 

0.3 

0.3 

0.3 

0.4 

0.3 

^Aa. 

0.763 

0.737 

0.718 

0.699 


0.743 


0.009 

0.010 

0.009 

0.008 


0.016 

^Ab. 

0.100 

0.119 

0.130 

0.145 


0.113 


0.009 

0.010 

0.009 

0.008 


0.016 

(b . 

0.137 

0.144 

0.152 

0.156 


0.144 


fixed 

fixed 

fixed 

fixed 


0.018 

XAa . 

0.74 

0.50 

0.48 

0.35 


0.50 


fixed 

fixed 

fixed 

fixed 


0.03 

^Ab. 

0.83 

0.60 

0.54 

0.42 


0.60 


fixed 

fixed 

fixed 

fixed 


0.03 

VAa. . 

0.33 

0.33 

0.33 

0.33 


0.33 


fixed 

fixed 

fixed 

fixed 


fixed 

VAh . 

0.41 

0.41 

0.41 

0.41 


0.41 


fixed 

fixed 

fixed 

fixed 


fixed 

o- (mag). 

0.0071 

0.0082 

0.0081 

0.0071 


0.0098 

^obs. 

732 

738 

739 

740 


6129 


^The second row for each parameter indicates the uncertainty (standard deviation), or 
whether the parameter was held fixed in the fit. All solutions adopt a mass ratio q = 0.7266 
from Table 10. 
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Table 10. Orbital solutions for V1061 Cyg. 


Parameter Spectroscopic only Combined 


Adjusted quantities from inner orbit (Aa and Ab) 

Pa (days). 2.3466487 zb 0.0000049 2.34665473 zb 0.00000035 

7 (kms-i). -7.84 ± 0.25 -5.67 ± 0.14 

ii'Aa (kms-i). 87.96 ± 0.34 87.83 ± 0.20 

J^Ab (kms-i). 120.54 ± 0.73 120.87 ± 0.70 

BA. 0 (fixed) 0 (fixed) 

(deg). 

Ta (HJD-2,400,000)'* . 52015.9050 ± 0.0018 52015.89295 ± 0.00058 

Adjusted quantities from outer orbit (Aa+Ab and B) 

^’ab (days). • ■ ■ 5786 ± 76 

ifA (kms-i). 5.37 ± 0.26 

eAB. • ■ ■ 0.469 ± 0.032 

iAB (deg). • ■ ■ 112.0 ± 4.9 

LJA (deg). • ■ ■ 97.2 ± 2.8 

Hab (deg). • ■ ■ 27 ± 44 

Tab (HJD-2,400,000)‘’ .... ■ ■ ■ 48545 ± 55 

Other adjusted quantities 

Ao* (mas). • ■ ■ +1.2 ± 1.9 

A5 (mas). +0.40 + 0.67 

A/i* (mas/yr). • ■ ■ -2.5 ± 3.5 

A/ij (mas/yr). • ■ ■ —4.0 ± 1.9 

Att (mas). • ■ ■ —0.26 ± 0.20 

Derived quantities 

Xb (kms-i). 12.85 + 0.53 

Ko-C (days). • ■ ■ 0.01458 + 0.00048 

/r* (mas/yr). +17.4 + 3.5 

(mas/yr). • ■ ■ +35.5 + 1.9 

TT (mas). • ■ ■ 5.99 + 0.20 

A^Aa (Mq) . (1.274 + 0.018)/sin3iA 1-282 + 0.015 

^Ab (Mq) . (0.9297 + 0.0096)/sin3 lA 0.9315 + 0.0068 

Mb (Mq). • ■ ■ 0.925 + 0.036 

q = MAb/MAa. 0.7297 + 0.0053 0.7266 + 0.0042 

Ma (Mq) . 2.213 + 0.021 

AfA + MB(M0). 3.139 + 0.046 

UA (Rq) . (9.666 + 0.038)/siniA 9.681 + 0.031 

“ab (ins'®). ■ ■ ■ 55.3 + 1.8 

flAB (AU). • ■ ■ 9.235 + 0.082 

“ph (™a'‘5). • ■ ■ 9.55 + 0.54 


^Time of primary eclipse. 
'^Time of periastron passage. 
















































Table 11. 


Physical parameters for V1061 Cyg. 


Parameter Primary Secondary Tertiary*^ 


Mass (M 0 ). 

... 1.282 

± 

0.015 

0.9315 

± 

0.0068 

0.925 

± 

0.036 

Radius (R©). 

... 1.615 

± 

0.017 

0.974 

± 

0.020 

0.870 

± 

0.087 

logs'. 

... 4.129 

± 

0.011 

4.430 

± 

0.018 

4.525 

± 

0.088 

Temperature (K). 

... 6180 

± 

100 

5300 

± 

150 

5670 

± 

150 

logL/Lo. 

... 0.533 

± 

0.030 

-0.173 

± 

0.052 

-0.153 

± 

0.066 

My (mag). 

... 3.456 

± 

0.074 

5.457 

± 

0.079 

5.23 

± 

0.16 

nsini (kms“^). 

36 

± 

2 

20 

± 

3 

2 

± 

3 

UsyncSini (kms“^). 

34.8 

± 

0.4 

21.0 

± 

0.4 




Distance (pc). 

... 166.9 

± 

5.6 

166.9 

± 

5.6 

166.9 

± 

5.6 

m — M (mag). 

... 6.113 

± 

0.073 

6.113 

± 

0.073 

6.113 

± 

0.073 


^The radius, \ogg, and luminosity of the tertiary were inferred from My, temper¬ 
ature, and bolometric corrections by Lejeune, Cuisinier & Buser (1998). 

















Table 12. Parameters for eclipsing systems in the mass range of VI061 Cyg. 



Mass 

Radius 

logs 

Tuff 

My 

Period nsini^ logLx^ 

Star 

(Mq) 

(R©) 

(cgs) 

(K) 

(mag) 

(days) (kms~^) (erg s~^) 


V1061 Cyg Aa. 

1.282 

± 

0.015 

V1061 Cyg Ab. 

0.9315 

± 

0.0068 

V1061 Cyg . 

0.925 

± 

0.036 

FL Lyr A. 

1.221 

± 

0.016 

FL Lyr B. 

0.960 

± 

0.012 

RW Lac A. 

0.928 

± 

0.006 

RW Lac B. 

0.870 

± 

0.004 

HS Aur A. 

0.900 

± 

0.019 

HS Aur B. 

0.879 

± 

0.017 

Sun^ . 

1 

.000 


1.615 

± 

0.017 

4.129 

± 

0.011 

6180 

± 

100 

0.974 

± 

0.020 

4.430 

± 

0.018 

5300 

± 

150 

0.870 

± 

0.087 

4.525 

± 

0.088 

5670 

± 

150 

1.282 

± 

0.028 

4.309 

± 

0.020 

6150 

± 

100 

0.962 

± 

0.028 

4.454 

± 

0.026 

5300 

± 

100 

1.188 

± 

0.004 

4.257 

± 

0.003 

5760 

± 

100 

0.949 

± 

0.004 

4.409 

± 

0.004 

5560 

± 

150 

1.004 

± 

0.024 

4.389 

± 

0.023 

5350 

± 

75 

0.873 

± 

0.024 

4.500 

± 

0.025 

5200 

± 

75 


1.000 4.438 5780 


3.456 ± 0.074 

5.457 ± 0.079 

5.23 ± 0.16 

2.3467 

2.3467 

36 ± 2 

20 ± 3 

2 ± 3 

30.12 

3.95 ± 0.09 

5.37 ± 0.10 

2.1782 

2.1782 

30 ± 2 

25 ± 2 

30.19 

4.48 ± 0.09 
5.12 ± 0.13 

10.3692 

10.3692 

2 ± 2 

0 ± 2 

< 29.6 

5.23 ± 0.08 

5.68 ± 0.08 

9.8154 

9.8154 


< 29.2 

4.83 


2 

26.4-27.7 


^Measured rotational velocities for the eclipsing systems are all consistent with synchronous rotation. Rotational velocity measurements for 
HS Aur are not available, but Popper et al. (1986) have shown that the lines are sharp. Synchronous velocities for the primary and secondary are 
5.2 kms“^ and 4.5 kms~^, respectively. 

'^The values listed are for each system as a whole. 

^Values for the radius and logg are inferred from other properties; see §5.3. 

"^The range in the X-ray luminosity of the Sun represents the change during the activity cycle (Peres et al. 2000). 
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Fig. 1.— Sample spectra of V1061Cyg (left panels) and corresponding cross-correlation 
functions (right) showing a broad-lined (Aa) and a sharp-lined (B) component. The Julian 
date, year, and orbital phase {ip) are indicated. The tertiary star (Ab) is not obvious in the 
1987 spectrum, but is in many of the more recent observations. 
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Fig. 2.— Radial velocity residuals for the primary and secondary of VlOGlCyg from a 
preliminary double-lined orbital solution, and radial velocities for the tertiary. The trends 
indicate the system is a hierarchical triple. 
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Fig. 3.— O — C diagram of all available times of eclipse for V1061 Cyg, from a preliminary 
linear ephemeris derived from a fit to the photoelectric (hlled circles), visual (open circles) 
and photographic measurements (triangles) for both the primary and secondary minima. 
Oscillations suggesting a third body are obvious, particularly in the recent data. 
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Fig. 4.— URSA U-band photometry for V1061Cyg along with our htted light curve, and 
residuals from that £t. 
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Fig. 5.— URSA U-band photometry for VlOblCyg around the primary minimum, along 
with our model light curve. Residuals from the £t are shown in the bottom panel. 
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Fig. 6.— URSA U-band photometry for V1061 Cyg around the secondary minimum, along 
with our model light curve. Residuals from the £t are shown in the bottom panel. Note that 
the vertical scale in the top panel is different from Figure 5. 
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Fig. 7.— BVRI observations for VlOblCyg from Gettysburg, along with our calculated 
light curves and residuals from those hts. The light curves in different passbands are shifted 
vertically for clarity. Patterns in the residuals near phase 0.45 may be due to spots (see 
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Fig. 8.— HIPPARCOS photometric observations for VlOOlCyg along with onr light curve 
£t using EBOP. The bottom panels show enlargements of the regions around the minima. 
This fit was used to derive the fractional luminosity of star B in the Hp passband (see text). 
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Phase 

Fig. 9.— Radial velocity observations for stars Aa (filled circles) and Ab (open circles) in the 
eclipsing binary, with our htted orbit. The motion of the binary in the wide orbit has been 
subtracted. Phase 0.0 corresponds to the primary eclipse, and the dotted line represents the 
center-of-mass velocity of the system. 
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Fig. 10.— Radial velocities for V1061 Cyg in the outer orbit, and 0 — C timing residuals as a 
function of time. The curves represent our best-fit model that uses velocities simultaneously 
with the eclipse timings and HIPPARCOS observations, (a) Velocities for the third star 
(hlled circles) and the center-of-mass of the inner binary (crosses). The latter is computed 
at each date from the weighted average of the two stars, after removing the motion in the 
inner orbit. The center-of-mass velocity is indicated with the dotted line, (b) 0 — C diagram 
of the more recent eclipse timings based on the linear ephemeris in Table 10. Open circles 
are for visual measurements, and hlled circles for photoelectric/CCD measurements. The 
older photographic timings have much larger error bars and are not shown, for clarity. The 
dates of the HIPPARCOS observations are indicated, and happen to cover periastron passage 
(arrows). 
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Fig. 11.— Path of the center of light of VlOhlCyg on the plane of the sky, along with 
the HIPPARCOS observations (abscissae residuals). See text for an explanation of the 
graphical representation of these one-dimensional measurements, (a) Total motion resulting 
from the combined effects of parallax, proper motion, and orbital motion according to the 
global solution described in the text. The arrow indicates the direction and magnitude of 
the annual proper motion, (b) Residual orbital motion after subtracting the parallactic and 
proper-motion contributions. The plus sign represents the center of mass of the triple, and 
the direction of motion in the orbit (retrograde) is indicated by the arrow. Periastron is 
shown with an open circle. 
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Fig. 12.— Comparison between the observed properties (error boxes) of the primary and 
secondary of V1061 Cyg (mass, radius, absolute visual magnitude) and model isochrones by 
Yi et al. (2001) and Demarque et ah (2004) for solar metallicity and three different ages, as 
labeled. 






55 




Mass (Mg) 


Fig. 13.— Age and metallicity range of theoretical models by Yi et al. (2001) and Demarqne 
et al. (2004) that fit the measnred properties of the primary and secondary components 
of VlOOlCyg. (a) Shaded areas represent age-metallicity combinations consistent with the 
measured mass, radius, and absolute visual magnitude of each star. The greyscale level is a 
measure of how well each model matches the exact values of M, R and My (darker areas 
represent a better £t); (b) The shaded area represents the projection onto the mass/radius 
diagram of all models indicated in the top panel that £t the primary star within its error 
boxes. None of the models come close to reproducing the secondary radius; (c) Same as 
above, for the mass/My plane. The secondary star is well fit by the same models that 
match the primary. 
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Fig. 14.— Observations of V1061 Cyg (error boxes) compared against model isochrones by 
Baraffe et al. (1998) for three different valnes of the mixing length parameter, as labeled. 
Solar metallicity is assnmed, and the age of each isochrone has been tnned to £t the properties 
of the primary as closely as possible within errors (2.6 Gyr for auL = 1-0, 3.1 Gyr for 
«ML = 1-5, and 3.4 Gyr for auL = 1-9)• The model with ctml = 1-0 (solid line) is seen to 
provide a mnch closer match to the radins of the secondary than the one with auL adjnsted 
to fit the Snn. 
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Fig. 15.— Properties of V1061 Cyg and three other eclipsing binaries having one star closely 
matching VlOhlCyg Ab in mass, compared against the Yonsei-Yale models by Yi et al. 
(2001) and Demarque et al. (2004). For each binary system the model isochrone best htting 
all observed properties is shown: Z = 0.0118 and age = 10.7 Gyr for RW Lac, Z = 0.024 
and age = 10.5 Gyr for HS Aur, Z = 0.026 and age = 2.4 Gyr for FL Lyr, and Z = Zq 
and age = 3.4 Gyr for V1061 Gyg. The secondary star in FL Lyr is signihcantly larger than 
predicted, as is V1061Gyg Ab (see text). The Sun (©) is also shown for reference, along 
with the third component of V1061 Gyg (§5.3). 
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Fig. 16.— Mass of the convective envelope for main-sequence stars as a function of stellar 
mass, based on evolution models by Siess, Forestini & Dougados (1997) and Siess, Dufour & 
Forestini (2000) for solar composition and the estimated age for V1061 Cyg of 3.4 Gyr. The 
envelope mass is shown in solar masses (top) and as a fraction of the stellar mass (bottom). 
According to these models stars of 0.4 Mq and lower are fully convective. In the inset the 
region around the masses of V1061 Cyg is shown on a larger scale, and the envelope mass of 
both stars is indicated. 




